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Abstract 

We consider particular modification of the free-field representation of the form factors in the 
Bullough-Dodd model. The two-particles minimal form factors are excluded from the construction. 
As a consequence, we obtain convenient representation for the multi-particle form factors, establish 
recurrence relations between them and study their properties. The proposed construction is used 
to obtain the free-field representation of the lightest particles form factors in the $1,2 perturbed 
minimal models. As a significant example we consider the Ising model in a magnetic field. We check 
that the results obtained in the framework of the proposed free-field representation are in agreement 
with the corresponding results obtained by solving the bootstrap equations. 

1. Introduction 

As it is known two-dimensional statistical models in their critical points are described by the so called 
minimal models of the conformal field theory pQ . Away from critical points the scaling region can be 
described by relevant perturbations of the fixed point action. The corresponding models can be referred 
to as the perturbed minimal models. The perturbations destroy long-range correlations of the critical 
model and the associated quantum field theories are usually massive. However in some cases an infinite 
number of integrals of motion survive. In particular, the minimal models perturbed by the one of the 
following primary operators, $1,2, $2,1, $1,3, are known to be integrable [2] [3]. Besides, in [4 h is 
shown that the $1,5 perturbation of non-unitary minimal models is also integrable. The $1^3 perturbed 
minimal models can be described as certain restrictions of the sine-Gordon model [SUB]. The other cases 
are related with particular quantum group restrictions of the Zhiber-Mikhailov-Shabat model [7J [HI H] - 

The off-shell behavior of these models can be studied in the framework of the form factor approach [S] . 
In particular, correlation functions can be reconstructed using the spectral density representation. The 
fast rate of convergence of the spectral series for all distance scale [10] allows one to calculate them quite 
accurately ignoring the multi-particle form factors. Besides, in the recent years another application 
of the form factor approach is deeply developing in the framework of studying a specific class of the 
two-dimensional non-integrable field theories. The non-integrable field theories can be treated as a 
particular perturbations of the integrable ones [TT] and mult i- particle form factors are significant in the 
study of these models. A particular example is provided by the Ising model whose off-critical behavior 
can be considered by using the two different integrable models. The recent progress in the studying 
of the magnetic deformation of the Ising model at non-critical temperature can be found in [12 . The 
perturbation theory near another integrable point is more complicated (13] . Some numerical data were 
obtained in [141 115) by means of the truncated conformal space approach. 

The mentioned relation between different integrable models provides an effective method of the form 
factors calculation in the perturbed minimal models. For example, the iS-matrix of the particles in the 
$1,3 perturbed minimal model originates from the breathers S'-matrix of the corresponding sine-Gordon 
model whose off-shell behavior was considered in a number of papers [HI [TT]. However for the $1,2, 
$2,1 and $1,5 perturbed minimal models the situation is quite different because of complexity of the 
S- matrix of the Zhiber-Mikhailov-Shabat model [7J[H]. In this paper we consider a very particular class 
of the form factors which can be obtained from the form factors of the Bullough-Dodd model [HI [19] [20] . 
The convenient method of the multi-particle form factor calculation is provided in the framework of the 
free-field representation [21] . This method was successfully applied for finding multi-particle form factors 
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in many integrable models [HJ [531 12H US . The free- field representation for the Bullough-Dodd model 
was proposed in [37J [25] . However in some cases particular modification of the free- field representation 
is more convenient. In this representation the two-particle minimal form-factors are excluded from the 
construction. In f2TTlf5U] this method was applied for finding form factors of descendant operators. In this 
paper we consider another its feature. Namely, this representation possesses simple analytical properties 
and can be used to obtain convenient free-field representation for the lightest particles form factors in 
the $1,2 perturbed minimal models. 

As it is known in the free-field representation the calculation of the multi-particle form factors reduces 
to a combinatoric procedure. The A-particle form factors in the Bullough-Dodd model and, consequently, 
in the corresponding perturbed minimal models are given by the sum of 3 N terms. In opposite to 
conventional construction in the proposed free-field representation Wick's averaging procedure generates 
functions with simple analytic properties. This feature allows us to obtain particular recurrence relations 
between the multi-particle form factors and prove the reflection properties for the form factors of the 
exponential operators that was suggested in |31j . Also, we prove that the form factors satisfy the quantum 
equation of motion. The certain restrictions of the imaginary coupling Bullough-Dodd model correspond 
to $1,2, $1,5 or $ 2j i perturbations of the M. PtP i minimal models p3[HlH]. ^ n ^is paper we consider the 
first possibility. The S-matrix of the lightest particles in these models 05| correspond the S-matrix of 
the Bullough-Dodd model being analytically continued to the imaginary values of the coupling constant. 
This allows us to propose the free-field representation for the lightest particles form factors in the $1,2 
perturbed minimal models. 

As a significant example we consider the Ising model in a magnetic field (IMMF). This model 
correspond to the $1,2 perturbation of the A^3,4 minimal model and is known to be related with 
the Eg algebra [3J. It is a challenging problem to obtain the free-field representations associated with 
this algebra directly. The remarkable connection between IMMF and the Bullough-Dodd model was 
established in [27J HH]- We present the free- field representation for the lightest particle of the IMMF 
explicitly and compare the calculations performed in the framework of the free-field representation with 
the corresponding results obtained in [T31 [33] by solving the bootstrap equations. We checked that these 
results are in full agreement. 

This paper is organized as follows. In section [3J we review the general features of two-dimensional 
integrable quantum field theories and specify the notation we use. In section[3Jwe present the modification 
of the free-field representation by excluding the two-particle minimal form factors from the construction. 
We discuss the features of the proposed representation and obtain particular recurrence relations between 
the form factors. By using this relations we prove the reflection properties of the form factors of the 
exponential operators and show that these form factors satisfy the quantum equation of motion. In 
section|3]we consider quantum group restrictions of the complex Bullough-Dodd model which correspond 
to the $1,2 perturbations of the A4 P , P ' minimal models and propose the free-field representation for the 
lightest particles form factors in these models. As an example of the proposed construction the Ising 
model in a magnetic field is considered in detail. 

2. Scattering theory 

Let us consider the relativistic scattering theory containing n sort of the particles A ai a = 1, . . . , n with 
masses m a . We introduce the notation A a {9) for the particle A a having the rapidity 9, 

p° a = m a cosh9 a , p 1 a = m a smh9 a , (2.1) 

where p° a , p\ are components of the two-dimensional momentum p a such that the mass-shell condition 
is satisfied, i.e. p\ = rr? a . The states 

|A ai (0i),...A a „(0„)) m(out) (2.2) 

form a basis of the asymptotic in(out)-states which is assumed to be complete. The rapidities of the 
states are supposed to be ordered, that is, 9\ > . . . > 9 n for incoming asymptotic states and 9\< ... <6 n 
for outgoing ones. 

The S'-matrix is defined as the unitary transformation connecting "in" and "out" asymptotic states. 
Let us assume that all the particles A a have different masses. This assumption, in particular, implies that 
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all the particles are neutral, A a — A a , where A a denotes the corresponding anti-particle. Let us restrict 
ourself to the case of the diagonal scattering only. Then the ^-matrix is diagonal in the basis (|2.2[) 

14,! (fil), • • • 4 n (0n))in = S ai ...an (9i, ■ ■ ■ , 9 n )\A ai . . . A a „ (0 n ))out, (2.3) 

An infinite number of the conserved currents in the integrable theories forces the scattering processes to 
be completely elastic, i.e. the final state contains the same number of particles with the same momenta as 
the initial one. Therefore, the n-particle S'-matrix element in (|2.3[) factorizes in the product of n(n — l)/2 
two-particle ones, 

S ai ... an (Oi,... ) n )='[[Sa i a j (Oi,O j ). (2.4) 

i<j 

The amplitudes S ab {9i,02) are meromorphic functions in the rapidity difference 9 — 612 — 9\ — 62, real 
at Re 9 = and satisfying the equations 

S ab (9) = S ab (nr - 9), S ab (9)S ab (-9) = 1, (2.5) 

which express the crossing symmetry and the unitary conditions of the theory respectively. From 
these conditions it follows that the amplitudes S a b{9) are 27ri-periodic functions, which are completely 
determined by positions of their zeros and poles in the "physical strip" < Im 9 < n. The poles on the 
strip are located at Re 9 = 0. Simple poles with positive residue correspond to "bound state" particles in 
the s-channel of the A a A b scattering, whereas those with negative one correspond to "bound states" in 
the w-channel. The basic bootstrap requirement is that the "bound states" must belong to the same set 
of the particles A a , a — 1, . . . , n. The simple poles with positive residue of the two-particle amplitudes 
S ab (9) at 9 — iu c ab represents the particle A c with mass 

m c = ml + ml + 2m a m b cos ul b , < h e(0,7r), (2.6) 

in the s-channel of the A a A b scattering. In the vicinity of the bound state pole the two-particle amplitude 
becomes 

S o6 (0~i„^)~i£|£-, (2.7) 

where Y c ah can be referred to as the on-shell three-particle coupling constant. 

The off-shell behavior of the integrable models can be studied in the form factor framework approach. 
For the local operator O(x) its form factors are defined as the matrix elements of this operator in the 
basis of the asymptotic states, 

-Ff;... a' m a„(#m) • • •>0'll#l) • - • > #n) = a' t . . .a' m {9\, ■ ■ ■ , 9' m \0 (0) 1 6\ , . . • , 6»„) Ql . . .a„ • (2.8) 

The crossing symmetry condition of the S'-matrix (|2.5p suggest that the form factors can be expressed 
in terms of the following matrix elements, 

F ° 1 ... an (9 1 ...9 n ) = (vac\O(0)\9 1 ...9 n ) ai ... an , (2.9) 

where (vac\ is the vacuum state (without particles) of the scattering theory. 

Let us review general conditions providing a systematic way to reconstruction of the form factors. For 
the scalar operator 0(x) the requirement for relativistic invariance implies that its form factors depend 
on the rapidity differences 9i — 9j only, while for the operator with the spin s one gets the condition 

+ A, ...,#„ + A) = e sK FZ.. an {0i, • • • . On)- (2.10) 

The form factors (12.9[) satisfy a set of requirements, called the form factor axioms, which are supplemented 
by the assumption of maximum analyticity [5]. They are the following, 

I. Watson's theorem, 

F °a k a k + 1 ...(9l, ■ . . ,9 k k +l, ■ ■ ■ 9 n ) = S(9 k+ i - 9k)F®a k+1 a h ...{Ql) • • • ) 9k+l,9 k , ■ ■ ■ n ). (2.11) 
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II. The crossing-symmetry condition, 

F°a 2 ...a n (0U^ ■ ■ -A) = Fa 2 ...a nai (02, ■ • -A, <?1 + 2™). (2.12) 

III. The kinematical pole condition, 

- * lim AO' - 6)FZ ax ... an {0' + in,9,9 1 ,..., 9 n ) = 

v — ytf 

n 

= (l-]JSaa i (6-8jj)FZ... an (9i,...,9 n ). (2.13) 

IV. The bound state pole condition, 

- % lim (9' - 9)F°^ an (9 1: ...,9' + iu c bd , 9 - tu b cd , ...,6 n ) = 

v — }o 

= ^LFZ.. d ...aA0 1 ,...,9,...,9 n ), (2.14) 

where u c ab = n - u c ab . 

These axioms do not specify the dependence of the form factors on the operator O(x) and further 
requirements are necessary to identify the form factors of a specific operator among all solutions to the 
form factor axioms. This is in general a non-trivial problem. However there are the two useful criteria. 
The first one argued in [33] restricts the asymptotic behavior of the form factor. For the scaling operator 
0{x) of the scaling dimension 2Aq the form factor growth for large values of rapidities is bounded by 

F°.. Mn (9 u ...,9 N )<e A °^, as|0(|->oo. (2.15) 

The second criterion essentially depends on the normalization of the form factors. For later use let us 
choose the following normalization, 

C. On (0i, • • • , On) = (O)fZ.. an (0u- • • , On), (2.16) 

where (O) is the vacuum expectation value of the local operator O. This choice of the normalization is 
consistent with the spectral representation of the two-point correlation function 

(O(x)O(0)) = £ £ 1^ . . . f^...^, . . . A)| 9 e-£"* COshej (2-17) 

™=0{a„} 

Further we refer to the function f® . . . , 0jv) defined in (|2.16p as to the non-normalized form 

factor. In these notation the second criterion proposed in |34j can be formulated in the following way. 
The form factors of the spineless exponential operators possess the factorization property, namely, 

fai...a n (0l + A, . . . , 6» m + A, 9 m+ i, . . .9 N ) = /oi„.o m (01> • • • >0m)f am+1 ... an {0m+l ■ ■ ■ 9 n ), (2.18) 

as A — > oo for any m 6 (0, . . . , n). This property had already been noticed to be satisfied in the solutions 
of many models and is believed to be a distinguishing property of exponential operators. 

The solution to the system of the equations (|2.11l) - (12.14l) can be represented in the convenient form 
by using minimal two-particle form factors. The minimal two-particle form factor, denoted by R a b{9) 1 
is a solution to the equations (|2.11[) and (|2.12[) being specified to the case n = 2, 

R a b(0) = S ab (9)R ba (-9), R ab {m + 9) = R ah {m-9). (2.19) 

In order to determine an unique solution to these equations one imposes further restrictions, namely, 
Rab{0) is an analytic function in the "physical strip" and has nor zeros nor poles in this range. By using 
these notation the n-particle solution to the form factor axioms can be represented in the following form, 

fZ...a n (0 1 ,...,9 n ) = jZ... an (e 6 \...,e s -) J] R aiaj (0i-0j). (2.20) 

l<2<j<n 
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Here the function j£ a (e Sl , . . . , e 9 ™) is a symmetric and 27ri-periodic function in the variables &i 
containing all expected kinematical and bound state poles. This function represents a solution to the 
second pair of the form factor axioms (I2.13[) and (|2.14p . Further we refer to these equation as to the 
bootstrap ones. The whole set of the form factors can be uniquely determined by the form factors which 
involve "fundamental particles" only. We referred to the particle Ai as to the fundamental one if all 
other particles of the model A a , a — 1 , . . . , n can be obtained as the "bound states" of some number of 
A\. As the form factors involving "fundamental particles" are calculated the other form factors can be 
obtained by using the residue equation (|2.14[) . Therefore we restrict our attention to the calculation of 
the multi-particle form factors. 

Further we consider form factors of the exponential operators only. In this case in framework of the 
free-field representation one can reduce their calculation to a combinatoric procedure |21j . In the next 
section we briefly review the free-field representation for the form factors of the exponential operators in 
the Bullough-Dodd model. 



3. The Bullough-Dodd model 

In this section we briefly describe a scattering theory of the Bullough-Dodd model [18] . This model is a 
two-dimensional integrable quantum field theory defined by the Euclidean action 

Sbd = J d 2 x (^{dvV? + K^ b<p + Ste"^)) , (3.1) 

where b is the coupling constant and [i is the regularized mass parameter which is related to the mass 
m of the single bosonic particle A in the spectrum of the model in the following, 

m 2V ^) ( ^ l + 2h2 ) \ & ( 2^(l + 6 2 /2) x sfc 

r(i + A)r(^)V r(-2*) ) I r(-B»/2) ) ' (3 ' 2) 

Here we introduce the convenient notation 

Q = b- 1 +b. (3.3) 

The Bullough-Dodd model possesses an infinite number of local conserved charges of the odd integer 
spin s except multipliers of 3, 

s = 1,5,7,11,13,.... (3.4) 

The integrability of the model implies that the n-particle S'-matrix factorizes into the n(n — l)/2 two- 
particle scattering amplitudes, 

_ tanh 1(0 + gf) tanh 1(0 - j|) tanh 1(0 - 
U tanh |(6I- 2hl) tanh 1(0+ |^) tanh \{9 + l ' J 

For real values of the coupling constant b the S'-matrix has a simple pole at = 2z7r/3 corresponding to 
the bound state represented by the particle A itself in the scattering processes 

AxA^A^AxA. (3.6) 

The space of local operators of the model consist of the exponential operators 

V a (x) = e av{x) (3.7) 

and their descendants. We are interested in the form factors of the exponential operators only. Let us 
label them in the following, 

F Q (0 1 ,...,0 Ar ) = (vac\V a {x)\e 1 ,...,e N ). (3.8) 

By using free-field representation one can obtain a set of functions which satisfy the form factor axioms 
together with the requirements (|2.15p and (j2.18[) . Consequently these functions can be identified with the 
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form factors of the exponential operators. We briefly review the free-field construction in the Bullough- 
Dodd model [35]. We do not need the detailed description in terms of the auxiliary free- field. Thus, let 
us formulate the result only. 

Consider a pair of the operators A(9) and A^ 1 (6) and define the braces ((. . .)) by 

«A(<?)»=1, 

({A°'(6')A*(9))) = [R(9-d')f' 7 , a',a = ±l, 

where R(9) is the two-point minimal form factor of the Bullough-Dodd model [THIHO]- It is convenient 
to use the following parametrization for this function, 

R(9)=M(b) fl 3 ' (3.10) 

Here we introduce the functions g a (9) defined by 

/ f°°dtcosh^Dl ( l7r -0) t x 

^ = cxp ( 2 / T^fc sin -^r-)' (3 - u) 

The normalization constant J\f(b) is chosen such that R(9) — > 1 as 9 — > 00, namely, 

/ P°° dt cosh I sinh A sinh -A^ \ 

N(b)=eU-A d ± 6 . , . 1 3Qb ■ (3.12) 

V J t sinh t cosh 7j / 

By construction the braces ((. . .)) with arbitrary number of the operators A(8) and A~ 1 (9) factorizes 
into the pair braces 



((A^(8 N )...A^(9 1 )))^ [] ((A ff *(^)A^(^)». (3.13) 

l<i<j<N 

With these notation Lukyanov's generators are given by 

T(9) = p ( 7 A(0 + I ) + 7 - 1 A- 1 (A - y ) + h : A(9 + ^)A(9 - y ) :), (3.14) 
where the ^-independent constants 7, /i, p are the following, 

7= -zexp(^(4V2a-6 + 6- 1 )), 

, n . 7r(fo - fo- 1 ) 

fe = 2sm 6Q ' (3-15) 



/ sin f exp r 2 / 

sin ^ «n ^ " 7 



00 dt cosh I sinh $ sinh 3^ 



3Q5 - * sinhtcoshf 

Notice that these notation differs from those one of [28 by the substitution b — > b^ 1 . The free-field 
construction possesses the duality transformation 6 <H- which is also a symmetry of the S- matrix. 
However, this substitution is applied to adopt the conventional notation of the conformal field theory. 
As a result, the form factors of the exponential operators can be represented as follows, 

F a {9 x , ...,0 N )= (e^rpu ...,9 N ) = (e a «>)((T(9 N ) . . . T(9 1 ))), (3.16) 

where ^-independent factor (e av ) is the vacuum expectation value of the exponential operator in the 
Bullough-Dodd model [3T]. By using (I3.9[) one can show that the non- normalized form factors f a (9\, . . . , 6n) 
in (|3.16[) possess the following parametrization, 

f a (0i,...,9 n )=p N J N , a (xi,...,x N ) Yl R(9 i -9 j ). (3.17) 

l<i<j<N 

Here Jn,h(xi, ■ ■ ■ , xn) is a symmetric rational function in the variables x\ = e 9i with proper kinematical 
and dynamical poles. Notice that in the proposed free-field representation the calculation of the form 
factors reduces to a combinatoric procedure. In the next section we obtain the free-field representation 
for the functions Jn,o.(xi, ■ ■ ■ >%n) alone. As a consequence we propose particular recurrence relations 
between them. 
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3.1. Free-field representation for the functions Jn,o,(xi, • ■ • ,x n ) 

We follow the guideline of Lukyanov's construction of the free- field representation [55] • Let us specify 
the differences. We consider the Heisenberg algebra generated by a countable set of the elements a*, 
n € Z, n and s = ±, 0, with the following commutation relations, 

r ^ ./i , • /(n — 1W& 7r mr\ . f(n— l)nb w\ . .„ . „. 

K,<J - -^KSg + 6 " t) sm (Sg^ + eH' '"' (3A8) 

and the matrix A a,(T is given by 




exp 



(f), (3.1.) 



where the raw s = +, — , and the column s' = +, — , 0. Let us define the exponential operators 

X s (z) = :expV^z™ : , s = ±,0. (3.20) 
^— ' n 

Here the normal ordering is not needed for the operators A (z) since the corresponding generators 
commute. Following we can determine the proper Lukyanov's generators by using the operators p. 201) 
and identify the exponential fields with the projectors on the corresponding Fock modules. It is convenient 
to formulate the result as a following prescription. The function Jn,o,{xi, ■ ■ ■ , %n) can be represented as 
a matrix element, 

Jn,o,(xi,...,x n ) = (t(xi) ...t(x N )), (3-21) 

where 

t(z) = 7A+(z) + j- x \-(z) + h\°{z) . (3.22) 
Here 7 and h are defined in (I3.15[) . In (I3.21j) we use the notation 

(...) = <0|...|0) (3.23) 

for the corresponding matrix elements and the vacuums are defined by the relations 

<|0)=0, (0|al„ = 0, s = ±,0, forn>0. (3.24) 

The calculation of the matrix elements multi-particle matrix elements is performed using Wick's averaging 
procedure with the following rules, 

A ± (a;)A ± (a;') = : A ± (a;)A ± (x') : , 
\°(x)\°(x>) = /(-^) : A (z)AV) : , 
A°(a;')A ± (a;) = A ± (x)A°(a;') = /(^w* 1 ) : A ± (a;)A (x') : , 
X-(x')X + (x) = X+(x)X-(x') =/(^ w )/(^ 2 ) :X+(x)X-(x'): , 

where the function f(x) is given by 

f(x) = — — : — • (3.26) 

x + x 1 — 1 

The free-field representation (|3.2ip is very similar to Lukyanov's representation (|3.16p . Its main advantage 
is that redundant factors of the two-particles minimal form factors are excluded from the construction. 
Therefore, calculation of the matrix elements involve Wick's averaging represented by rational functions 
in the variables X4. Notice that the averaging procedure p. 251) in the ./V-particle matrix elements generates 
3 N terms. The particular terms may possess poles at the points x% = XjUj ±2 corresponding to bound 
state poles, at Xi = —Xj corresponding to kinematical poles, and redundant ones at ±1 . The 



(3.25) 
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presence of the redundant poles complicates calculation of the form factors. However, in Appendix [Al we 
prove that their contributions in the whole expression for the matrix elements vanish. 

The simple analytic structure of the matrix elements allows us to obtain recurrence relations between 
them. Consider the function Jjv+i, a (^, xi, ■ ■ ■ >%n) a $ an analytic function in the variable z depending 
on the parameters x±, . . . , xjv- One can separate the contribution of the poles from the regular part. 
Since the residues of this function can be evaluated explicitly (see Appendix [5| the only unknown is the 
regular part. However, for the exponential operators the regular part of this function can be determined 
using the factorization property (|2.18[) . We present details of the construction in Appendix [Bj As a 
result we can formulate the following proposition, 

Proposition 1. The recurrence relations 

JN+l t a{z,%l, ■ ■ ■ ,%n) — Jl,aJN,a(xi, . . . ,X N ) + 

N 

+ \^ — — K n (xi,. . . ,X N )J N -i Jxi, . . . , jf„, . . . , X N ) + 



n— 1 
N o 



-B+(X!, . . . ,X N )J N , a (xi, . . .,X n LO, . . .,X N )- 



(3.27) 



Z — X n U! 
n—1 



N 

- V xB~(x\, . . . , X N ) JjV,a(a;i, . . . , XnU!' 1 , . . . ,X N ). 

— X n 0J l 

n—1 

together with the initial condition 

J ha = J = 4sin(^^)cos(^(2y2a - b + b" 1 )) , (3.28) 

uniquely define the set of functions Jn,o,{xi, ■ ■ ■ , xn) of exponential operators. The functions K n (x\, . . . , Xjv) 
and (xi, . . . , xn) are given by 



h{h 2 - 1 

LU — US 



B*(x 1 ,...,x N )= =T" u ) 

L) — UJ 1 ± 1 \ Xi J 



i^n 



As an example of possible applications of the recurrence relations (I3.27[) in Appendix[C]we prove that 
the form factors of the exponential operators satisfy the quantum equation of motion. Another example 
is provided by the so called reflection relations |31) . 



3.2. Reflection relations for exponential operators 

The Bullough-Dodd model can be interpreted as the perturbed Liouville theory in the two different 
ways, with either the first exponent in (|3.1[) or the second one is taken as the perturbing operator. Let 
us introduce the notation 

Q L = ^- b + V2b, Q' L = ^ + ^ 2 (3-30) 

There is a conjecture based on the features of the Liouville theory that the form factors of the exponential 
operators in the Bullough-Dodd model satisfy the relations [3Tj . 

{vac\e a *\e 1} ...0 N ) = B,(a)(vac\e^ L - a ^\0 u . . . , 9 N ), 
(vacle-^ie^ ...,9 N ) = i?'(a)(t>ac|e ( ~ Q ' L+a) ^i, ...,9 N ), 

where the explicit expression for the reflection amplitude R(a) is presented in |31j and R'(a) is obtained 
from that by the substitutions /i — > 2fi b — > b/2. The vacuum expectation values possess the same 
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relations. Consequently, for the non-normalized form factors defined in (|2.16[) the reflection relations can 
be presented in the following form, 

f a {6 1 ,...,e N ) = f^- a {e 1 ,...,e N ), 
r a (e 1 ,...,9 N ) = r^ +a (e 1 ,...,9 N ). 

By construction f|3 . 1 7[) the non-normalized form factors only depend on the parameter a through the 
functions Jn,o,{xi, ■ ■ ■ , %n) whose dependence on this parameter is completely determined by the initial 
conditions to the recurrence relations (|3.28[) . One can easily check that the functions J\ >a indeed satisfy 
the reflection relations. This proves (|3.32[) . 

3.3. An explicit calculation of the form factors 

From the representation (13.21)) it follows that the functions Jjsf, a (xi, . . . , xjy) are symmetric and 2iri- 
periodic functions in the variables Xi = e 9i . Consequently, these functions can be expressed in the 
convenient form by means of symmetric functions. A basis in the space of the symmetric functions in N 
variables is provided by the elementary symmetric polynomials <jn{x\, . . . , a; at) defined by the generating 
function 

N N 

Y\_(x + = yix N ~ % <Ti(xi, . . .,x N ). (3.33) 

i=l i=0 

It is convenient to separate the poles contribution prescribed by the recurrence relation (|3.27p into an 
overall factor. Then the functions Jn,o.(xi, ■ ■ ■ > %n) can be parametrized as 

J N , a {xi,...,x N ) = A N (a 1} ...,a N ) TT - ■ rrr~, ; — ST' ( 3 - 34 ) 

l<i<]<N ( Xl + X ^ X i + XiX i + X V 

where Ajv(o"i, . . . , cat) are symmetric polynomials in the variables ai{x\, . . . , xn)- Taking into account (12. 10)) 
one can show that for spinless operators the total degree of this polynomial is equal to 3N(N — l)/2 
while the partial degree in the each variable Oi is bounded by the the condition ()2 . 1 5[) and could not 
exceed 3(N — 1). In the framework of the free-field representation p.21[) we compute the polynomials 
Ajv(<ti, . . .,a N ) up to N = 4, 

Ai = J, 

A 2 (eri,cr 2 ) = J 2 a\- J (J + h-h 3 )a 2l 
A 3 (a u a 2 ,a 3 ) = J(J 2 + (-3 + h 2 )(-l + h 2 ) 2 + J(h - fi 3 ))^ 
+ J(J 2 + (-3 + h 2 )(-l + h 2 ) 2 + J(h - h 3 ))<j 1( j 2 

- Jh{J + h- /i 3 )(4 - bh 2 + h 4 )aia 2 al (3.35) 

- J 2 (J + h- h 3 )axa 2 

- J 2 {J + h-h 3 )<jt<j 2 a 3 

+ J(-J 2 + 3(-l + h 2 ) 2 + Jh(4: - 5h 2 + fe 4 ))of of ff 3 
J <J 1 a 2 , 



and explicit expression for the four-particle symmetric polynomial A4(<7i, a 2 , a 3 , a 4) is presented in 
Appendix [DJ 

As a simple check of the proposed expressions let us consider form factors of the exponential operators 
which appear in the action (|3.1D . The general expression of the stress-energy tensors trace ©(x) 
which is compatible with the quantum equation of motion is given by the linear combination of these 
operators [35]. Taking into account the stress-energy tensor conservation law one can show that the 
A-particle form factors of ©(x) has to be proportional to the combination of symmetric polynomials 
Gi&N-i f° r N > 2. Using explicit expressions of the three- and the four-particle form factors we checked 
that the form factors of the corresponding exponential operators indeed possess this property. 
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4. The $i )2 perturbed minimal models 



(21 

Consider an imaginary coupling afBne Toda theory based on the twisted affine Kac-Moody algebra A 2 , 
which is obtained from (|3.1I) by the substitutions b — ¥ i/3, /i — > —fx, 



S, 



cBD — 



(^(^) 2 -Me^-2e-^)). (4-1) 



This model can be referred to as the Zhiber-Mikhailov-Shabat or the imaginary coupling Bullough-Dodd 
model. It is not unitary theory since its Hamiltonian is not hermitian. However in [7] it was shown that 
the model is reducible for certain values of its coupling constant. Let us briefly describe the restriction 
of the complex Bullough-Dodd model. By using the action (|4.ip it is possible to determine non-local 
conserved charges that commute with the Hamiltonian. These charges generate the quantum affine 
algebra U q (A^') with 

q = e^ 1 (4.2) 
and fix the ^-matrix of the model up to a scalar factor El IB] - The restriction proceeds in the following 

(2) 

way. The quantum algebra U q {A 2 ) contains two subalgebras U q (sl2) and t/ g 4(sZ 2 ) and both of them can 
be used to obtain quantum group restrictions of (|4.ip . Consider the first possibility. Let the parameter 
2/3 2 takes the following rational values 

2/3 2 = - (4.3) 

where p and p' are relatively primary integers such that p' > p > 1. In this case the Hilbert space 
of the imaginary coupling Bullough-Dodd model can be consistently truncated to the representations of 
Uq{sh)- The details of the construction can be found in [7J[5]. The restricted theory is known to coincide 
with the $1,2 perturbations of minimal models M PtP ' [JJ whose formal action can be represented in the 
following form, 

Syu , +A [ d 2 x4> h2 (x). (4.4) 



The relation between the corresponding coupling constants /i and A can be established using the results 

of ng, 

4V r 2 (i-2/? 2 )r(4/? 2 ) 

(4/3 2 - l) 2 r 2 (2/3 2 )r(l-4/3 2 )' 1 ' 

Let us briefly introduce conventional notation used in the conformal field theory. The central charge of 
the minimal model A4 P . P > is given by 

c=l + 6Q 2 , (4.6) 

where Ql was introduced in (I3.30[) and the substitution b = i/3 is assumed. The primary operator content 
of the minimal model ftA P:P > is given by a set of the degenerate fields $ m .„, with m = 1, ... ,p — 1 and 
n = 1, . . . ,p' — 1, whose conformal dimensions can be parametrized in the following, 

A m ,n — &m,n(Q L ^m,n); i^-^) 

where 

(n- 1) ,-.„ (to- 1) 1 

(4.8) 

The normalization of the primary operators is fixed by the two-point function, 

($ m ,„(x)$ m ,„(0)) = |x|- 4A — \ (4.9) 

The primary operators of the perturbed minimal models relates with the specific exponential fields of 
the imaginary coupling Bullough-Dodd model 015]. Namely, the following exponential operators, 

V ai Jx) =e a ^ x \ (4.10) 

commute with particular generators of the C/ g (sZ 2 ) and can be identified with the primary operators $i jn 
of the perturbed minimal models. 
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The fundamental particle of the complex Bullough-Dodd model is a kink triplet. The relation between 
masses of the lightest kinks and the coupling constant A is established in (33]. Under the action U q (sl2) 
the fundamental kinks form either singlet bound states generating a series of breathers or triplet ones 
corresponding to higher kinks. After analytical continuation to the imaginary values of the coupling 
constant and quantum group restriction the S-matrix of the Bullough-Dodd model can be treated as 
the lightest particles S-matrix in the $1,2 perturbed minimal models [32 . It develops two additional 
poles in the "physical strip" corresponding to heavier particles. We propose an identification between 
the form factors of the Bullough-Dodd model and the lightest particle form factors of the corresponding 
minimal models. One can check that after analytical continuation to the imaginary values of the coupling 
constant the form factors (|3 . 1 Tf) satisfy the form factor axioms. In particular, the two-point minimal form 
factor develops two additional poles in the "physical strip" whose positions correspond to the poles of 
the S-matrix of the lightest particles and can be explicitly determined by using the following functional 
relations, 

9i+ a (0) = 9- a (0), 

ia -cosh (4-11) 



la {9)g^ a {9)=V a {9) = 



2 cos 2 2f 



Further, we restrict our attention to the non-normalized form factors only. In the free-field representation 
the correct normalization is fixed by the vacuum expectation value of the corresponding operator. 
However, perturbed minimal models possess several vacuums and calculation of their contribution to the 
vacuum expectation values is an open problem. As a result let us propose the following representation 
for the lightest particles form factors of the primary operators in the $1^2 perturbed minimal models 

f^l(e 1 ,...,e N ) = r-(9 1 ,...,e N ) . (4.12) 

2b 2 =— p/p' 

Notice that the form factor in the right hand side of these relation admits the free-field representation (|3.21l) 
Being analytically continued to the imaginary values of the coupling constant the minimal two-particle 
form factors develop additional poles. Therefore, the poles emerge in Lukyanov's representation (|3.16[) 
due to Wick's averaging procedure (|3.9[) . However, in the proposed free-field representation (|3.21l) these 
redundant factors are excluded from the construction. The averaging procedure generates functions with 
the same analytic structure as in the Bullough-Dodd model. In the next section we consider the Ising 
model in a magnetic field as an significant example of the proposed free field representation. 

4.1. Ising model in a magnetic field 

The model (|4.4[) with (p,p') — (3,4) describes the Ising model at critical temperature in non-zero 
magnetic field (IMMF). The minimal model ^3,4 is characterized by its central charge c = 1/2 and 
contains two non-trivial primary operators $1.2 and $1,3. Their conformal dimensions are given by (|4.9p . 
namely, Ai^ = 1/16 and A13 = 1/2. The Ising model can be defined in a different way, namely, as a 
scaling limit of the Ising lattice model at its critical point. Therefore, one can define the spin density 
operator a(x) and the energy density operator e(x) as the scaling limit of the corresponding lattice 
operators. Assuming the normalization (|4.9[) for these operators one can identify them with primary 
operators of the minimal model, namely, <j(x) = $>i^{x) and e(x) — $1,3(2;). The coupling constant A 
defined in (|4.5[) can be referred to as external magnetic field or, more precisely, as its scaling limit. 

The spectrum of the model consist of eight different species of self-conjugated particles A i: i = 1, . . . , 8. 
Under analytical continuation to imaginary value of the coupling constant, 2b 2 = —3/4, the S'-matrix of 
the Bullough-Dodd becomes the S'-matrix of the lightest particles in the IMMF, 

tanh 1(0 + gf) tanh \{9 + 2 -f) tanh \{9 + jf ) 
1U ' tanh±(0- tanh *f )tanh±(0- § )' K ' ' 

The additional poles that emerges in the "physical strip" after analytical continuation of the S-matrix 
correspond to heavier particles. We denote the particle corresponding to the pole at = 27ri/3 by A 2 
while those one corresponding to the pole at = iir/15 by A3. As shown in [3] the full set of the 
two-particle amplitudes S o b(0), where a, b = 1, ... 8, can be reconstructed starting from the two-particle 
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amplitude (I4.13[) . Besides the bootstrap structure of the model fixes masses of all particles, 

mi =m, mo = zmcos-, 7713 = 2m cos — , 

' 5 30 

77T 27T 7T 

m 4 = 2m 2 cos— , m 5 = 2m 2 cos — , m 6 = 2m 2 cos — , (4.14) 

oU to ou 

it 7n f 2ir 

7717 = 47772 COS — COS , J718 = 4m 2 COS — COS — , 

5 30 5 15 

Notice that under analytical continuation of the Bullough-Dodd model and quantum group restriction 
not all integral of motion survive. Instead of Q3.4p the IMMF possesses integral of motion with 

,s = l,7,ll,13,17,19,23,29, mod 30. (4.15) 

The integers (|4. 1 5[) are are exactly the exponents of the Lie algebra Eg, repeated modulo 30. In addition 
the number of particles is exactly the rank of Eg. 

4.2. Form factors in the IMMF model 

The form factors of the IMMF can be calculated in the framework of the proposed free-field representation. 
It is instructive to consider the analytic continuation of the form factors (|3.17[) in detail. The analitical 
continuation of the two-particle minimal form factor can be performed by using relations (|4.1ip . 

where we introduce the notation 

Rii(0)=rf(iyjf) So(% 2 / 3 (%2/5(%i/i5(#)- (4-17) 

One can show that this function satisfy the form factor axioms (|2.11l) . (|2.12l) and has nor zeros nor 
poles in the "physical strip" . Therefore we refer to this function as to the minimal two-particle form 
factor of the lightest particles A\ in the IMMF. The poles representing by the functions VaiO) in (|4.16|) 
correspond to the A 2 and A 3 bound state poles of the two-particle amplitude (|4. 13[) . 

The analytic continuation of the function Jn,o{ x i, ■ ■ ■ ,%n) is straightforward. By construction this 
function satisfy the form factors axioms (|2.14l) and (|2.13|) at the dynamical poles, 9 = 27ri/3, and 
kinematical ones respectively, has a correct large rapidity behavior (|2.15[) and possesses the factorization 
property (|2.18p . In the parametrization A3. 17)) the lightest particle form factors in the IMMF can be 
represented in the following form, 

/-(*„.. .,.„>- „»^„ fa ,. ..,„, || „ 18) 

where the substitutions 2b" — —3/4 is assumed. Notice that the functions Jjv,ai n {x±, . . . , xjy) in (|4.18[) 
are given by the corresponding matrix elements (I3.21j) . Therefore, one can obtain multi-particle form 
factors of the lightest particles either by using Wick's averaging procedure or by means of the proposed 
recurrence relations. We checked that both these methods up to four-particle form factors reproduces 
the results of [13J . where the form factors were obtained by solving the bootstrap equations. Due to 
bootstrap structure of the model the form factors involving different species of particles can be obtained 
from (|4.18p by means of the residue equation (|2.14l) . For example, in the proposed representation the 
computation of the four-particle form factor f t \'™ 2 ( a; i> ^3) x i) reduces to the averaging procedure 
of the corresponding matrix element J§. ai n (xi, x%, x±, 25CJ 3 / 5 , 25W~ 3 / 5 ). Notice that in this model the 
correct normalization of the form factors of primary operators can be obtained by using the vacuum 
expectation values of the corresponding exponential operators in the Bullough-Dodd model |31) , 

($1,2) = -1.27758 . . . A 1/15 , ($ 1)3 ) = 2.00314 . . . A 8/15 . (4.19) 

These values are in good agreement with numerical calculations |37j and have been verified for the 
operator $13 with a higher precision in [Tl]. As a result we show that the proposed construction of the 
free-field representation can be used for the mult i- particle form factors computation in the IMMF. 
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5. Conclusion 



Let us summarize the results. In the framework of the free-field representation we construct a space of 
solutions to the form factor axioms in the Bullough-Dodd model. The proposed free-field representation 
differs from the conventional by excluding the minimal two-particle form factors from the construction. 
As a consequence, the calculation of the multi-particle form factors reduces to the computation of 
particular matrix elements. The iV-particle matrix element is given by the sum of 3 N terms whose 
calculation can be performed by using Wick's averaging procedure generating rational functions in the 
variables Xi — e 9i only. The particular terms in the sum may possess non-physical poles. However, we 
consider the operator algebra in detail and prove that the contribution of these poles in the whole sum 
vanishes. 

A simple analytic structure of the proposed matrix elements allows us to obtain explicit recurrence 
relations between them. By using these relations we prove that the form factors of the exponential 
operators are consistent with the quantum equation of motion and satisfy reflection properties. Besides, 
in the framework of the free-field representation we compute the form factors of the exponential operators 
up to four-particle one. 

We consider the quantum group restrictions of the imaginary coupling Bullough-Dodd model. In 
contrast to Lukyanov's construction of the free-field representation where additional poles necessarily 
arises from the two-particle minimal form factors contribution, the proposed construction can be easily 
used to obtain the free-field representation for the lightest particle form factors in the $1.2 perturbed 
minimal models. Since we exclude the minimal form factors from the construction, the computation of 
the matrix elements is performed by Wick's averaging procedure generating the rational functions with 
the same set of poles as in the Bullough-Dodd model. 

A particular example is provided by the Ising model in a magnetic field. It is a challenging problem to 
obtain the free-field representations associated with the Eg algebra directly. The remarkable connection 
between this model and the Bullough-Dodd model was established in [27J HH] • I n this paper we present 
the free-field representation for the lightest-particle form factors in the IMMF in convenient form. As 
a consistency check we verify that the results obtained in the framework of the free-field representation 
are in full agreement with those ones obtained by solving the bootstrap equations. 
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Appendix A. Analytic properties of the function Jn,o,( x 1) ■ ■ ■ , x n) 

Consider the function Jjv+i,a(z, x\, ■ ■ ■ , %n), which is a symmetric function in the variables z,x±, . . . , Xff. 
Therefore it is sufficient to determine its analytic properties with respect to the variable z alone considering 
the other ones as a parameter. Due to the form factor axioms (|2.13[) and (|2.14p one can expect that the 
only poles of the function Jjv+i,a(z, x\, . . . , xm) are those at the points 

z = Xiio ±2 , z = —Xi, (A.l) 

where the former poles correspond to the rapidity difference 9 = 2ni/3 and, therefore, to the bound state 
poles, while the later one gives 9 = in and corresponds to the kinematical pole. Let us prove that the 
matrix element 

(t(z)t(x 1 )...t(x N )) (A.2) 

being considered as a function in the variable z has simple poles at (|A.1[) as well. Besides we determine 
the corresponding residues explicitly. We consider the operator product expansion 

t(z)t(x). (A.3) 
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By using (|3.25p one can chow that its simple poles are those at the points 



z = x^ 1 , z = xuj ±2 , z = —x. 



(A.4) 



Let us consider these poles consequentially. 



I. The poles at z = xu> . As it follows from (IA.1I) these poles correspond nor dynamical nor 
kinematical ones. Let us prove that the corresponding residues of the matrix element (|A.2|) vanish. 
The singular part of the operator product expansion (|A.3[) is given by 



(z - xu ±l )t{z)t{x) = ±xuj ±1 rh 2 (: A°(xw ±1 )A°(a;) : - : A T (xuj ±1 )X ± (x) :) 

where we introduce the following notation, 

r ( Z \ +1 tl 2 — 1 

Res / I — I = ±xuj r, r = -. 

z=xui ±1 \X/ UJ — UJ 1 



The operators that appears in the right hand side of (|A.5|) satisfy the following relations 

(: A°(a;cj ±1 )A (x) : X a {x')) = (: \ t (xlu ±1 )X ± (x) : X a {x')), a = ±,0. 
Therefore, from (|A.5[) together with (|A.7|) we obtain 

(z - a;„u; ±:L )(£(z)t(xi) . ..t(x n ) . . .t(x N )) 



(A.5) 



(A.6) 



(A.7) 



(A. 



II. The dynamical poles at z = xoj ±2 . At these poles the singular part of the operator product 
expansion is given by 



- XL0 ±2 )t(z)t(x) 



±xto ±2 hrt(xLJ ±L ) 



(A.9) 



where we introduce the operator t(x), 
t(x) = 7 : X°(xuj +1 )X + {xuj^ 1 ) : : A~ {xlo +1 )X°(xuj' 1 ) : +h : X~ {xuj +1 )X+ {xuj- 1 ) : . (A.10) 



Notice the similarity between this operator and those one defined in (|3.22[) . In every matrix element 
containing the operator t(z) one can replace it by t(z) taking into account that 



N 



n=l 



(t(z)t( Xl )...t(x N )) = l[f^) (t(z)t(x n ) . ..t(x N )). 
Therefore from (|A.9[) and (jA.111) we obtain 
(z - x n u ±2 )(t(z)t(xi) . . . t(x n ) . . . t(x N )) 



(A.11) 



Z — X n UJ A 



where 



= x n u ±2 B± { Xl ,..., xnMx^ 1 )^) . . .tlxrf. . . , t(x N )), (A.12) 
B±(xi, ...,x N ) = hr]l f^^ 1 ) • (A.13) 



III. The kinematical pole at z — —x. The singular part of the operator product expansion at this pole 
is the following 

(z + x)t(z)t{x) =xf(io 2 )r(s(xio- 3 / 2 )~s(xuj^ 2 )) (A.14) 

z— — x V / 



where the operator s(x) is given by 

s(x) =: X+(xlu-^ 2 )X-(xlo 3 / 2 ) 



(A.15) 
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Being inserted into the matrix element this operator simply produces an overall factor, 

N 



( S (z)t( Xl )...t(x N )) = J{f(-u 1 ' 2 )f(-u- 1 ! 2 '){ti{x l )...t{x N )) (A.16) 
Consequently, from (|A.14[) and (|A.16I) we get 
(z + x n ){t{z)t{xi) . . . t{x n ) . . . t(x N )) = 

Z — — X 

= x n K n (xi, . . .,x N )(t(x!) . . .tix^j . . .,t(x N )), (A. 17) 

where 

*•<*• *»> -f^im^-')^-') - n/(t^)0))- (A.i8) 

Thus we prove that the matrix element (|A.2I) as a function in the variable z has simple poles at (|A.1|) 
and the corresponding residues (|A.12|) and (|A.17|) can be represented as a matrix elements of the less 
number but the same operators t( X i). 

Appendix B. The regular part of the function Jn,o,( x 1i ■ ■ ■ ,^iv) 

Consider the function Jn,o,{z, X) as an analytic function of the variable z depending on the parameter 
X = {xi, . . . , a; at}. It is possible to separate the contribution of the poles from the regular part, 

J N+l , a {z,X) = Jjx°+ 1<a {z,X) 
N 

+ J2^^K n (X)J N _ 1 (X n ) 

Z ~r X n 

n—1 

N 2 (B 1) 

+ V ^ B+(X)J Nta (x n LJ,X n ) 
^ z - X n CJ z 
n—1 

N _ 2 

- V* yB~(X)J N a (x n u~ 1 ,X n ). 

n—1 

where we introduce the notation X n = {xi, . . . , yt n , . . . , xn}. The function a {z, X) is a regular 

function in the variable z everywhere except the points z = and z — oo. Since the sum over poles is 
of order 0(z~ 1 ) as z -> oo, the asymptotic behavior of the function Jjv+i.a(zj X) as a function in the 
variable z is governed by J < ^ 1 a (z,X). The expansion (|B.1|) is applicable in the vicinity of zero. 
One can consider the expansion of the function Jn+i,o,(zi X) in the vicinity of the point z = oo, 



J N+1>a (z,X) = J ( °' +la (z,X) 
N _i 
- E il ^ K n (X)J N ^(X n ) 



n=l ~ 1 " n 
N -12 

E i n W -i ^ B+{X)J N>a {x n u +x ;X n ) 
/ , — — — ~B~ (X)JN :a (x n ui~ ; X n ), 

n=l " 



(B.2) 



where the function a (z, X) is regular everywhere except the points z — 0, oo. At the vicinity of the 

(0) 

N+l 



point z = the behavior of the function Jjv+i l0 («, X) is governed by jjv+i a (z, X 
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Notice that the expansions (|B.1[) and (|B.1[) are applicable for any kind of operators under consideration, 
i.e. not only for exponential but for descendant operators as well. At this point we restrict our attention 
to the exponential operators only. As discussed in section [5] the form factors of the exponential operators 
should satisfy the factorization property. Consequently the only unknown functions J^+i a {z,X) and 

J^ +la (z,X) in the recurrence relations can be fixed easily. Indeed, consider the asymptotic of the 
function Jjv+i. a (z,X) as z — > and z — > oo. Using the factorization property (|2.18[) we get 

J^ + \jz,X) = J$l lia (z,X) = J ha J N , a {X), (B.3) 
where the one-particle form factor is given by 

Ji, a = 7 + 7~ 1 + h. (B.4) 

Taking into account explicit expressions for h and 7 we obtain the proposition [TJ 

In the remain part of this section we obtain the useful identity for the functions Jn,o,( x 1i ■ • ■ ,Xn)- 
Let us define the function 

N N 

D n , a (X) = Y,Kn{X)J N -i{X n ) + J2x n u 2 B+(X)J N , a {x n oj +1 ;X n )- 

n—l n—1 

N 

-Y J ^- 2 B-{X)J N AxnOJ- 1 :X n ). (B.5) 



n=l 



From the expansions (jB.ll) and (|B.2|) it follows that 

•4°!i,a(*>*) - = D N , a (X) (B.6) 

Taking into account (|B.3p we get 



D N , a = 0. (B.7) 

This equality provides a non trivial identity for the functions Jjv,o(^i) ■ ■ ■ ,%n) and we use it in Appendix [Cl 
to prove that the form factors of the exponential operators are compatible with the quantum equation 
of motion. 

Appendix C. The equation of motion for form factors 

In this section we prove that form factors in the Bullough-Dodd model are consistent with the equation 
of motion, 

ddip = 8V2TT^b(e V2bv -e'^) . (Cl) 



The derivatives of a field produce multiplication of its form factors by the components of the momentum 
according to the usual rule P M — > id^. Let us introduce the notation 

S n (X)=xl + ... + x n N , S n (z,X) = z n + S n (X), (C.2) 

where X — {x%, . . . , xn}- Let z — e e and x n — e n . The components of the momentum are given by 

777 77? 
P z (9,d 1 ,...,6 N ) = --S 1 (z,X), P M (9 ) 9 1 ,...,6 N ) = -S- 1 (z,X). (C.3) 



Therefore, 

^ Sl (z,X)S^(z,X)j a 
Let 

J , N+1 (z,X) = ±J N+1 , a (z,X) _ n . (C.5) 



vac\dd<p\6,6 1 ,...,e N ) = -—S 1 {z,X)S- 1 (z,X)—f a (9,6 1 ,...,d N ) . (C.4) 

a— 



a=0 



Using these notation we can represent the equation of motion in the following form 

S 1 (z,X)S- 1 (z,X)J' N+1 (z,X) - A(j^ +liVl6 - J N+h _ b/y #), (C.6) 



16 



where the constant A is given by 

A = - " v "7" (e v ^), (C.7) 



and we take into account that the vacuum expectation values of the corresponding exponential operators 
coincide [3T] 

( e v^) = ( e -WV5) = 1^! J_ (C.8) 

A* 48y3Q&sin(^)sin(3^) 

Let us prove the equation of motion (|C.6I) recursively. Suppose that this equation is valid for some value 
N. Taking derivatives of both sides of the recurrence relation (|3.27[) we get 

N 

S 1 (z,X)S- 1 (z,X)J' N+1 {z,X) = Y,x n K n (X n )S- 1 (X n )J' N _ 1 {X n )+ 

71=1 

N 

+ E ^^K n {X)S 1 {X n )S- 1 (X n )J' N _ 1 {X n ) + 
z -\- x n 

n—l 

N 

+ ^ x n u 2 B+(X)S-i(x n LJ,X n )J' N+1 (x n uj,X n )+ 

n=l 
N 2 

+ Y"] (X)Sx(x n oj, X n )S-i(x n u;, X n )J' N (xuj, X n )- 

^— ' z - x n uj z 

n—l 

N 

- ^ x n ^~ 2 B~(X)S-i(x n u~ 1 ,X n )J' N+1 (x n oj~ 1 ,X n )- 

n=l 

N _ 2 

~ E ^B-(X)S 1 (x n u- 1 ,X n )S-i(x n u- 1 ,X n )J' N (xu- 1 ,X n ). 



(C.9) 



n=l 



Here in the right hand side one can apply the induction hypothesis (IC6|) . Besides we use the following 
identities, 

N N 2 

E ^r^-K n {X)J N ^ a (X n ) + J2 - , Xn " > - B+{X)J Nta {x n u,X n ) 

N _ 2 

~ E <r - B-{X)J N , a {x n u-\X n ) 

= -J N+ x, a {-Sx{X),X) + Ji, a J N , a (X), (CIO) 

and 

N N 2 

V -^-K n {X)J N - lia {X n ) + V XnU 2 B+J N , a (x n u,, X n ) 

' V -X- T— ^ ' r — T„f,? Z 



, 2 + X r , 
n—l 



N _ 2 
X n Ul 



~ E 2B-(X)J N , a (x n u;-\X n ) 

n—l 

= J N+ l, a {z, X) - J ha J N , a (X). (C.ll) 

where we use (|B.6|) . Therefore (|C.9I) gets the following form 

S 1 (z,X)S-i(z,X)J' N+1 (z,X) = A{J N+lV 2 b (z,X) - J N+lib /^2(z,X)- 

-JN+i,V2 b (Si(X),X) + J N+lb/ ^{-S 1 (X),X)). (C.12) 
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Compare this equation with those of (|C.6I) . To prove the equation of motion for N + 1 it is sufficient to 
prove that the last two terms vanishes, 

- J N+ i,V2 b (Si(X),X) + Jn+w^i-StiX), X) = 0. (C.13) 

In fact, these terms vanish separately independently on a. To prove (jC.13l) notice that the function 
JN+i.a{—S\{X), X) is symmetric in the variables X4. Therefore, any of these variables can be chosen for 
z, i.e. 

Jn+\m{—Si,X) = Jn+i,o,{— Si(xn + i, Xj), xn+i, Xj). (C14) 

The left hand side is ^-independent, while the right hand side is Xj -independent. This proves that the 
function <//v+i,a(— S±(X), X) is constant in all variables X. Let us show that this function is equal to 
zero. Consider the limit xn — > 00. From the recurrence relation we obtain 

jv ^ 

JN+i,a{— Si(X); X) = Ji, a Jjv,q(-X") + - K n (X)Jff-i ia (X n ) 

n=l bl{X n ) 



N 



,2 



+ E g /"% B+(X)J N , a (x n u +1 ;X n ) (C.15) 
n=1 bi{x n u,X n ) 

N _ 2 

~ E c 5 B-{X)J N A^ 1 ;X n ) 

n=1 bi(x n uj \ X n ) 

Since the left hand side is a constant, we can calculate it in the limit xn — > 00. In this limit the only 
non-vanishing terms in the sums in the right hand side are those with n — N. Using the factorization 
property one can show that in the leading order in xm the function Jjv+i,a( — Si(X); X) — > 0. Since this 
function is a constant we prove (|C13|) . Consequently we prove that the form factors Jn,o,( x i> ■ ■ ■ >%n) 
indeed satisfy the equation of motion of the Bulloug-Dodd model. 

Appendix D. The four-point symmetric polynomial A 4 (x 1 , . . . , x 4 ) 

In this appendix we present coefficients of the polynomial A^eri, 02, 03, 04). Let us choose the following 
parametrization, 

18 

A 4 (0i,0 2 ,0304) = Ci j! fc,i0i02 cr 3 cr 4- (D- 1 ) 

i,j,k,l=0 
i+2j+3fc+4i=18 
i+j+k+l<9 

In this notation the coefficients Cij.k.i are the following 

c ,o,2,3 = - J 4 - J(-l + ft 2 ) (3 J - 24ft + 7J 2 h + AJh 2 + 7 Ah 3 - 6J 2 h 3 

- 13Jft 4 - 85ft 5 + J 2 ft 5 + 7Jh 6 + A5h 7 - Jh s - llh 9 + h 11 ) 

Co, 1,0, 4 = 

c ,i,4,i = Jh{A - bh 2 + h 4 )(-3 + J 2 + Jh + 7h 2 - Jh 3 - bh 4 + h 6 ) 
c ,2,2,2 = 2J 4 + J(-l + ft 2 ) 2 (-6J + 21ft - 5 J 2 h - hJh 2 - 

- 46ft 3 + J 2 h 3 + 6Jh 4 + 3Ah 5 - Jh 6 - 10ft 7 + h 9 ) 
Co, 3, 0,3 = 

c ,3,4,o = J 4 — J 2 {-1 + ft 2 )(-3 + Jh + Ah 2 - h 4 ) 

c ,4,2,i = — J 4 + J(-l + ft 2 )(-3J - 3ft + 2J 2 h + bJh 2 + 7h 3 - 2Jh 4 - 5h 5 + h 7 ) 
Co, 5, 0,2 = 

ci,o,3,2 = 3J 4 + J(-l + ft 2 )(9 J - 2Ah + 13 J 2 h + AJh 2 + 7Ah 3 - \2J 2 h 3 - 25Jh 4 

- 85ft 5 + 2 J 2 ft 5 + 14Jft 6 + 45ft 7 - 2Jft 8 - lift 9 + ft 11 ) 
ci,i,i,3 — J 4 + Jh(-1 + h 2 ){-27 + 7 J 2 - 40Jft + 33ft 2 - 6J 2 ft 2 

+ 74Jft 3 + 46ft 4 + J 2 h 4 - AAJh 5 ~ 93ft 6 + HJft 7 + 52ft 8 - Jh 9 - 12ft 10 + ft 12 ) 
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ci,i,5,o = - J 2 KJ + h- ft 3 )(4 - 5ft 2 + ft 4 ) 

ci,2,3,i = -J i + + ft 2 ) 2 (3J + 15ft + J 2 ft + 17 Jh 2 - 19ft 3 - 9Jft 4 + 4ft 5 + Jh 6 ) 
ci,3,i,2 = — 2 J 4 - Jft(-1 + ft 2 ) 2 (27 - 5J 2 + 5Jft - 42ft 2 + J 2 h 2 - Jh 3 + 17ft 4 - 2ft 6 ) 
ci,4,3,o = -J 4 + J 3 ft(-1 + ft 2 ) 
ci,5,i,i = J 4 - J 2 h(-l + ft 2 )(2J + ft - ft 3 ) 

c 2 ,o,o,4 = - J 4 - J(-l + ft 2 )(3J - 24ft + 7J 2 ft + 4Jft 2 + 74ft 3 - 6J 2 ft 3 

- 13 Jft 4 - 85ft 5 + J 2 h 5 + 7 Jft 6 + 45ft 7 - Jft 8 - lift 9 + ft 11 ) 
C2,o,4,i = - 3J 4 - J 2 (-l + ft 2 ) 2 (-9 - 5Jft - 5ft 2 + Jft 3 + 6ft 4 - ft 6 ) 

C2,i,2,2 = - 3J 4 - Jft(-1 + ft 2 )(-54 + 21 J 2 - 36Jft + 138ft 2 - 14J 2 ft 2 + 69Jft 3 

- 118ft 4 + 2J 2 ft 4 - 43Jft 5 + 38ft 6 + llJft 7 - 4ft 8 - Jft 9 ) 
c 2 ,2,o,3 = 2 J 4 + J(-l + ft 2 ) 2 (-6J + 21ft - 5J 2 ft - 5Jft 2 - 

- 46ft 3 + J 2 ft 3 + 6Jft 4 + 34ft 5 - Jft 6 - 10ft 7 + ft 9 ) 
c 2 , 2 ,4,o = - J 4 + J 2 (-l + ft 2 )(-3 - 4Jft + 3ft 2 + Jft 3 ) 

c 2 ,3,2,i = 3J 4 - J(-l + ft 2 )(-6J - 6ft - 3J 2 ft + Jft 2 + 13ft 3 + 2J 2 ft 3 + 7Jft 4 - 8ft 5 - 2 Jft 6 + ft 7 ) 
c 2 A0,2 = - J 4 + J(-l + ft 2 )(-3J - 3ft + 2 J 2 h + bJh 2 + 7ft 3 - 2 Jft 4 - 5ft 5 + ft 7 ) 

C 2 ,5,2,0 = 

c 3 ,o,i,3 = 3J 4 + J(-l + ft 2 )(9J - 24ft + 13J 2 ft + 4Jft 2 + 74ft 3 - 12 J 2 h 3 - 25Jft 4 

- 85ft 5 + 2J 2 ft 5 + 14Jft 6 + 45ft 7 - 2 Jft 8 - lift 9 + ft 11 ) 
C3A5.0 - J 4 - J 2 M + ft 2 )(-3 + Jft + 4ft 2 - ft 4 ) 

cs,i,3,i = 3J 4 + Jft(-1 + ft 2 )(-3 + 13 J 2 + 4Jft + 7ft 2 - 8J 2 ft 2 - 5Jft 3 - 5ft 4 + J 2 ft 4 + Jft 5 + ft 6 ) 
cs,2,i,2 = - J 4 + J(-l + ft 2 ) 2 (3J + 15ft + J 2 ft + 17Jft 2 - 19ft 3 - 9Jft 4 + 4ft 5 + Jft 6 ) 

t4 

C3,3,3,0 — J 

C3,4,l,l = - J 4 + J 3 ft(-1 + ft 2 ) 

c 4 ,o,2,2 = - 3 J 4 - J 2 (-l + ft 2 ) 2 (-9 - 5 Jft - 5ft 2 + Jft 3 + 6ft 4 - ft 6 ) 
C4,i,o,3 = Jft(4 - 5ft 2 + ft 4 )(-3 + J 2 + Jft + 7ft 2 - Jft 3 - 5ft 4 + ft 6 ) 

C4,l,4,0 = - J 4 + J 3 ft(-1 + ft 2 ) 

C4,2,2,i = - J 4 + J 2 (-l + ft 2 )(-3-4Jft + 3ft 2 + Jft 3 ) 

C4,3,0,2 - J 4 - J 2 (-l + ft 2 )(-3 + Jft + 4ft 2 - ft 4 ) 

cb,o,3,i = J 4 - ^ 2 (-l + ft 2 )(-3 + Jft + 4ft 2 - ft 4 ) 
C6,i,i,2 = - J 2 H J + ft - ft 3 )(4 - 5ft 2 + ft 4 ). 
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